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ABSTRACT 

Tidal disruption flares are differentiated into two classes - those which are 
sub-Eddington and those which radiate near the Eddington limit. Flares from 
black holes above ~ 2 x IO^Mq will generally not radiate above the Eddington 
limit. For a Schwarzschild black hole, the maximum bolometric luminosity of 
a tidal disruption is ~ LEdd(5 x IO^Mq), substantially below the Eddington 
luminosities of the most massive disrupting black holes 2 X lO^M©). 

Bolometric corrections to the spectra of the brightest flares are found to be 
large ~ 7.5 mag. Nevertheless, the brightest flares are likely to have absolute 
magnitudes in excess of -19 in V and -21 in U (in the absence of reddening). 
Because the spectra are so blue, K-corrections may actually brighten the flares 
in optical bands. If such flares are as frequent as believed, they may soon 
be detected in low or high redshift supernovae searches. The He H ionizing 
radiation produced in the flares may dominate that which is produced by all 
other sources in the centers of quiescent galaxies, creating a steady state, highly 
ionized, fossil nebula with an extent of ~ 1 kpc which may be observable in 
recombination lines. 


Subject headings: accretion, accretion disks - black hole physics - Galaxy: 
center - galaxies: nuclei - quasars: general 


1. Introduction 


Many lines of indirect evidence suggest that massive black holes reside in the 
centers of a substantial fraction of galaxies (e.g. |Kees 1997| ; [Kormendy fc Kichstone 
1995|; and references therein). For example, remarkable stellar-dynamical studies of M32 
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( |van der Marel et al. 1997|) and maser studies of NGC 4258 (|Watson fc Wallin 1994 


pyiiyoshi et al. 1995|) show extreme central mass concentrations. Such concentrations are 
difficult to understand without a black hole. The most popular alternative explanation of 
extremely dense clusters, runs into problems at such high densities (e.g. [Goodman fc Lee 


|1989|) , although dynamical “mixtures” of a black hole with a small, dense cluster cannot 
be excluded (e.g. [van der Marel et al. 19971) . A nagging theoretical problem with black 
holes in quiescent galactic centers was that the black hole luminosities were too low for the 
expected levels of accretion; however the low luminosities can now be understood in context 
of low-efficiency, advection dominated models (e.g. [Narayan, Yi, fc Mahadevan 19951) . 


If black holes do, in fact, reside in the centers of galaxies, one clinching sign would be 
the tidal disruption of stars by the black hole (e.g. [Rees 198'5D . Importantly, the detection 
of a single tidal disruption event could even allow for a crude determination of the black 
hole mass (e.g. Loeb and Ulmer 1997). Such estimates would yield the black hole mass 
apart from any contribution from a dense stellar cluster. 


Tidal disruption events are rare (~ 10“^ per year per L* galaxy), bright, short (months 
to years), and seemingly unavoidable consequences of 10® — IO^Mq black holes in galactic 
centers. The general picture for tidal disruption of a star has been developed over the last 
20 years (e.g. Hills 1975; Young, Shields, & Wheeler 1977; Rees 1988). The picture has 
evolved from one in which the disruption rate was so high as to fuel quasars to one with a 
more modest rate in which disruptions are rare transients which may slowly grow a lower 
mass, ~ IO^Mq, black hole (Goodman & Lee 1989). 


If a star is scattered onto a near radial orbit (inside its loss cone) and at pericenter 
passes within the tidal radius, Rt, of the black hole, then the star will be disrupted as has 
been illustrated by a variety of numerical simulations (e.g. jlNolthenius fc Katz 1982| , [Evans 
fc Kochanek 1989] , [Laguna et al. 199^ , [Diener et al. 1997[ ). For the Sun, 


Rt — Rq{M]^yi /~ 25RsM^ 
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( 1 ) 


where Rs is the Schwarzschild radius, and Mg is the black hole mass in units of 10^Mq. 


As the star is ripped apart by the tidal forces of the black hole, the debris is thrown 
onto high eccentricity orbits with a large range of energy/period, where 




GMbhT?* 




( 2 ) 


and Rp is the pericenter of the star’s orbit ( [Lacy, Townes, fc Hollenbach 198^) . The 
distribution of mass as a function of energy is nearly constant as is borne out by simulations 
(e.g. [Evans &: Kochanek 1989| , [Laguna et al. 199'3[) . 
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Half of the debris is unbound with velocities around 5 x 10^ km/s. This unbound 
debris may possibly produce some long-term observable effects in the surrounding medium. 
Nucleosynthesis due to the extreme compression of a star may occur in a small fraction 
of tidal disruptions ([Luminet fc Pichon 1989| ; [Luminet fc Barbury 199U| ) and lead to 
observable enrichment within the Galactic center. In AGN and Seyferts, a highly-enhanced 
rate of tidal disruptions may produce enough unbound debris to form broad line and 
narrow line quasar-like regions (|Roos 1992|) . Khokhlov and Melia (1996) show that 
interactions between the high-velocity unbound debris and the ISM which may result in 
a supernovae-like remnant with properties similar to those of Sgr A East. In extremely 
low-angular-momentum encounters, the material could deposit over 10®^ ergs into the ISM, 
creating a very large supernova remnant. More typically, the energy deposited would be 
~ lO^o ergs. 

The bound half may create a bright flare as it accretes onto the black hole. The most 
bound material hrst returns to pericenter after a time. 
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( 3 ) 




The bound material returns to pericenter at a rate ([Rees 1988 , hhinney 1989|) given by 

1 M* / t \ 
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The peak return rate according to simulations by Evans & Kochanek (1989) is 


Peak Rate ~ 1.4 ( 
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which occurs at t 1.5G 


( 4 ) 


( 5 ) 


The scalings given above are for Rp Rs. When Rp ~ i?s, the tidal interactions 
Frolov et ah 1994|). These effects of the Schwarzschild metric can 


are stronger (e.g. 

be approximated in the pseudo-Newtonian potential of Paczyhski and Wiita (1980), 

4/(R) = —GM/{R — Rs). For instance, close to a black hole, the tidal radius, 

Rt,s = -^©(Afbh/MQ)^/^ -|- Rs, differs from that of the Newtonian potential. The deepest, 
non-captured orbit has Rp = 2Rg, so as is well-known, the maximum mass black hole that 
can disrupt a solar-type star is ~ 1 x lO^M©. As discussed in § IV, encounters more distant 
than the tidal radius may still strip signihcant matter from the star allowing disruptions 
for black holes up to ~ 2 x 10 ®Mq in the Schwarzschild case. Also altered is the range of 
energy of the debris (Eq. |^), 

GMbhR* 
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The increase in binding requires that Eqs. be modihed (for the pseudo-Newtonian 
potential, is replaced by (i?p — Rs)), and has the effect of reducing the minimum return 
time and of increasing the peak debris return rate. 


Another relativistic phenomenon, orbital precession, may lead to the crossing of debris 
orbits and may form strong shocks ([Rees 1988|) . The explicit evolution of the debris is 
complex and is dependent on the black hole’s properties as well as those of the star and the 
stellar orbit. As discussed by Kochanek (1994), at intersections between the debris streams, 
energy and angular momentum may be transfered. The computational complexities of the 
problem prohibit a direct numerical solution, but as an approximation, a large fraction of 
the debris can be said to circularize in a time. 


Air ^orbA 


(7) 


where riorb is a number greater than 1, and probably between 2 and 10. 


After, and perhaps while, the debris circularizes, it may accrete onto the black hole. 
The dynamical time at the tidal radius is short relative to the return time, so it is likely 
that the disk accretes quickly, resulting in a luminous flare. Cannizzo, Goodman, & Lee 
(1990) investigated the accretion of a tidally disrupted star in the thin disk case. 


Ulmer (1997) has investigated accretion in a thick disk case. The emergent spectra 
may have a strong dependence on the structure and orientation of the disk, as well as on 
the possible existence of winds. A starting point for models of the emergent radiation is the 
temperature associated with an Eddington luminosity emanating from the tidal radius 
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or from a small factor times the Schwarzschild radius: 
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(9) 


In either case, the temperature scaling with black hole mass is slight. Therefore, in the 
simplest models, all disruption flares would have comparable characteristic temperatures. 


A tidal disruption event has yet to be dehnitively detected, but there are a number of 
suggestive observations. First, a UV flare from the center of an elliptical was observed on 
the timescale of a year by Renzini et ah (1995), but this flare is 10“'^ as luminous as might 
be expected. Furthermore, the spectrum flattens in the near UV (~ 2000A) which would 
not be expected except perhaps in the case of a very strong wind (see § V). The possibility 
remains, as they suggest, that they have observed the disruption of only the outer layer 
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of a star. Second, Peterson and Ferland (1986) observed the optical brightening of an 
optically variable Seyfert along with the appearance of a strong, abnormally broad He II 
recombination line. The event is intriguing, but without a light curve or other dehnitive 
evidence, the event is not conclusive. As pointed out by Rees (1997), the most convincing 
disruption would be one which comes from a quiescent galaxy. Third, observations of 
variable double-peaked Hq, lines in spectrum of Seyfert NGC 1097 may be modeled as 
emission lines from the appearance of an elliptical disk (possibly formed through tidal 
disruption) close to a black hole, but other explanations exist (Eracleous et ah 1995; 
Storchi-Bergmann et al. 1995). 

The outline of the paper is as follows. In § II, a crude, but useful phenomenology of 
accretion disks is discussed. In § III, time scales in the tidal disruption event are presented. 
In § IV, these time scales are used to delineate regions of Eddington-level tidal disruption 
flares from the generally, less luminous sub-Eddington flares. In § V, characteristics of the 
emergent spectra from these flares are investigated. Conclusions are given in § VI. 


2. Accretion Disk Phenomenology 


In this section, the known relationships between luminosity and mass accretion rate 
are examined for a variety of accretion regimes which are important in determining the 
luminosity of tidal disruption flares. The mass accretion rate is given in terms of the 
Eddington accretion rate. Me = pEdd/ec^, where e is the efficiency of the accretion process 
and is taken to be 10%. 


Figure shows a schematic of luminosity versus mass accretion rate. The brightest 
known accretion disks have high mass accretion rates and are optically and geometrically 
thick, but thick disks are never much brighter than the Eddington luminosity, because at 
very high accretion rates, pressure gradients push the inner edge of the disk closer to the 
marginally bound orbit, so that less binding energy is released during accretion and more 
is carried into the black hole (e.g. [Paczynski fc Wiita 19801) . The luminosity increases 
logarithmically with mass accretion rate for thick disks ( [Paczynski 1980|) . Thin disks are 
less luminous, with luminosity scaling linearly with mass accretion rate (e.g. |Shakura~S 
Sunyaev 1978]). Slim accretion disks (jAbramowicz et al. 1988| jSzuszkiewicz, iVIalkan, fc 


Abramowicz 1995]) , may form a transition between the thin and thick accretion regime. 

At very low accretion rates, optically thin advection dominated models may exist and 
are even less luminous (Narayan & Yi 1995). They exist at a maximum accretion rate of 
~ o^Me, and have L ~ M^, so that the luminosity falls off faster with decreasing accretion 
rate. A possibility exists that at accretion rates very much higher than the Eddington rate 
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(and probably beyond those produced by stellar disruptions), the accretion may become 
more spherically symmetric and highly optically thick, so the radiative diffusion time could 
exceed the the infall time. In this case, much of the binding energy could be carried into the 
black hole (e.g. [Shvartsman 1971 ; [Ipser fc Price 1977| ; hlammang 1984 ) and the luminosity 
could actually fall below the Eddington limit. An additional complication is the possibility 
of beaming in the thick disk regime, but for simplicity and because of the absence of a 
preferred beaming model, beaming is set aside for the time being. 


The simplified view of accretion disks is that the luminosity increases with mass 
accretion rate until the Eddington mass-accretion rate, at which point the luminosity 
reaches the Eddington limit and remains at that level for all higher mass accretion rates. 
Accretion events with sub-Eddington luminosities occur in thin disks, and accretion 
proceeds in thick disks for all super-Eddington accretion rates. For sub-Eddington accretion 
rates, the luminosity falls off at least as fast as linearly with mass accretion rate. 


3. Derivation of Timescales 


Although many of the results derived in this section are already known, I briefly 
rederive them here for completeness and to scale all relevant parameters. The timescales 
considered are (1) the radiation time, the time to radiate, at the Eddington limit, the 
energy of the bound debris; (2) the circularization time, the time for a large fraction of the 
debris to circularize; and (3) the accretion time, the time to accrete a large fraction of the 
disk. Throughout the paper, m is defined as the mass of the star in units of the solar mass, 
r as the radius of the star in units of the solar radius. Mg as the black hole mass in units of 
10®Mq, as the Schwarzschild radius, 2GMbh/c^, and as the pericenter of the orbit in 
Schwarzschild coordinates. 


The main timescales of the problem can now be written in terms of the system 
parameters. 

The time to accrete all of the bound debris (which is half of the stellar mass) if the 
black hole radiates at the Eddington luminosity is 


O.SmMQC^e , 

Aad =- j ~ 21mMg eo.iyears 


( 10 ) 


for a hydrogen mass fraction, X = 0.74, and efficiency, e = 0.1. 

As discussed in §1, the debris are given a range of periods, and the circularization time 
is closely related to the the orbital time of the most tightly bound debris, fmin, (see Eq. ^ 
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and 1^: 


f ■ 5^ 

^cir ~ 


/ R \ ^ 

norb^min ~ O.llriorb ^^ears, 


( 11 ) 


where n is the small number of orbits necessary for circularization to occur. 

The time in which a large fraction of the matter in a disk can accrete can be written 




tKep(2i?p) 

3 X 10-^ 


aK^ 


^^^years, 


( 12 ) 

(13) 

(14) 


where h is the ratio of disk height to radius and is approximately one for thick disks, and 
txep is fhe rotation period in a Keplerian potential. In § ^ the geometrical thickness of the 
disk is discussed, but unless noted otherwise, it is the thick disk case that is of interest. 


4. Comparison of Timescales 

By considering the main timescales of tidal disruption events, the gross characteristics 
of the flares are investigated. Of particular interest is the luminosity and whether the flare 
is super-Eddington or sub-Eddington. The argument that will be made stems from the 
system’s time scales: trad, the time to radiate, at the Eddington limit, the energy of the 
bound debris; tcir, the time for a large fraction of the debris to circularize; and face the time 
to accrete a large fraction of the disk. 

Consider the three possibilities: (1) the accretion rate is so high that the material is 
drained from the disk nearly as fast as it circularizes, (2) the accretion rate is slower than 
circularization but faster than the radiation time so the material forms a reservoir which 
accretes in a thick disk, (3) the material accretes at a rate slower than it circularizes, and on 
a time scale longer than the radiation time, so the disk is not thick and the debris accretes 
in the less luminous form of a thin disk. 

Let us hrst consider the these scenarios for a solar-type star with pericenter at the 
tidal radius which is disrupted by a 10 ®Mq black hole. The parameter, riorb, related to 
the efficiency of circularization, is taken as 2. The only unspecihed parameter is a. For a 
IO^Mq black hole: (1) if a 1.5 x 10“^, then face ^ tch and a thick disk will form and drain 
as the material circularizes; (2) if 1.5 x 10“® ^ ^ 1-5 x 10“^, then ^acc ^ ^rad and 

a a reservoir of debris will form and accrete in a thick disk; and (3)if a ^ 1.5 x 10“®, then 





ta.cc ^ ^rad and a leservoir of debris will form and accrete in a thin disk. Therefore, nnless a 
is extremely small, mnch smaller than (~ 10“^ — 10“^) as is estimated from, for example, 
DN and FU ontbursts ( |(Jannizzo 1992| ; [Lin fc Fapaloizon 199fi| ; and references therein) and 
from nnmerical experiments of the Balbns-Hawley mechanism ( fiawley, Gammie, fc Balbns 
19951; [Stone, et ah 199'6|) , thick disks are the norm. Therefore, it appears likely that for a 


IO^Mq black hole, the tidal disrnption of a star prodnces a thick disk which will radiate 
near the Eddington limit (althongh it is still nncertain in what energy band the energy will 
be released). 


In hgnre these timescales are shown as a fnnction of black hole mass, again for the 
case of a solar-type star with a pericenter eqnal to the tidal radins. Becanse of the different 
scalings with black hole mass, the circnlarization time becomes longer than the radiation 
time at Mbh > 2 x lO'^eQ^iMQ. This Ending is independent of a as long as a is larger than 
~ 1.5 X 10“^. When the peak accretion rate falls below the Eddington rate above this 
critical mass, the accretion shonld occnr in a slim disk ( [Abramowicz et ah 1988|) or a thin 
disk rather than a thick disk. Althongh little has been said yet abont the spectra, it is clear 
that there may be observational difierences between the two classes. 


A similar conclnsion may be reached by considering when the peak retnrn rate of the 
stellar debris (eqs. |^) is snper-Eddington. For an accretion efficiency of 0.1, the maximnm 
mass is 1.5 x IO^Mq. This estimate is slightly smaller than that obtained by considering 
the timescales, becanse the circnlarization time is an nnderestimate, since a fraction of the 
matter retnrns in times longer than the circnlarization time. Nevertheless, both estimates 
point to a transition mass aronnd IO^Mq. For more massive black holes, the retnrn rate 
is lower, and the Inminosity, which for slim to thin disks scales ronghly linearly with mass 
accretion rate, will be lower than Eddington. 


The treatment above has be carried ont for a solar type star with pericenter eqnal to 
the tidal radins. Let ns first consider the resnlts for a less massive, main seqnence star. 
These stars will generally be denser with r ^ (e.g. [Kippenhahn fc Weigert 1990[) . 

Then, the transition mass will be ~ 2 X 


The last parameter that will be considered is the pericenter of the orbit which for 
disrnptions may be small as 2i?s in Schwarzschild coordinates. The npper radins, or the 
disrnption radins, i?d, is the radins at which a star in a parabolic orbit is barely disrnpted 
by the tidal field of a black hole. The disrnption radins is of order the tidal radins, and is 
given by 

Ra = qRt, (15) 

where the parameter, q, is between abont 1 and 3 and hides a nnmber of complexities 
inclnding relativistic corrections to the potential. Thronghont the paper, q is taken as 
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2. Young, Shields, and Wheeler (1977); Luminet and Carter (1986); Khokhlov, Novikov, 
and Pethick (1993); Frolov et ah (1994); and Diener et ah (1997) contain discussions 
of this parameter with a range of analytical approximations and numerical methods. 

The demarcation of is hazy as stars still may lose a fraction of their mass at larger 
pericenteric orbits. 

According to the theory of loss cones ([Frank fc Rees 1976t [Lightman fc Shapiro 19771) 
there is a critical orbital radius outside of which a star’s orbit diffuses in one orbit at an 
angle greater than the loss-cone angle. If stars which diffuse onto the loss cone primarily 
come from inside this critical radius, there will be an excess of events with i?p ~ R^. If 
stars diffuse primarily from large radii, where the diffusion angle is larger than the loss 
cone, one can use the “nan approximation” which Hills (1975) used to calculate the tidal 
disruption rate to show that the distribution of lost stars is constant with uj (gravitational 
focusing cancels the gain in cross section at larger uj). The true distribution of lost stars will 
be a combination of the two extremes mentioned, with an excess of events near i?p = 
and a tail towards R^ = 2Rs, corresponding to “plunge-orbits” in a Schwarzschild metric 
( [Novikov fc 'Fhorne 1973]) . The plunge orbits, in which the star is directly swallowed by 
the black hole, will not produce any tidal debris or accompanying flare. It is well-known 
that above ~ 2 x 10 ®Mq, the plunge radius surpasses the disruption radius, so no tidal 
disruption events can be produced. 

Tidal disruption events with smaller pericenters will have higher peak luminosities 
than those with Rp ~ Rt (see Eq. |^). Therefore, the actual transition mass discussed above 
is better regarded as a gradual rather than a sharp transition. In hgure the fraction of 
tidal disruption events with peak accretion rates greater than the Eddington-limit is shown 
for the case of a uniform distribution of Rp and stars with a Salpeter mass function between 
O. 3 M 0 and IMq. More realistic distributions of Rp, formed by integrating over a cusp, 
would provide a somewhat sharper transition. 

An interesting result of these scaling arguments is that if the luminosity is roughly 
limited at the Eddington luminosity, then there is a maximum luminosity that a tidal 
disruption event may produce. This luminosity occurs when the peak mass return rate 
is the Eddington rate and the pericenter is the minimum (2i?g) for a tidal disruption. 

The maximum luminosity corresponds to the Eddington limit for a black hole mass 
~ 5 X e^i Mq, i.e. 7 x 10^® ergs s“^ and 2 x IO^^Lq. The maximum luminosity is 
probably somewhat lower, because the estimate assumes that all of the gravitational 
binding energy of the returning debris is released as it returns to pericenter, whereas the 
circularization process is probably slower. Such extreme events would be rare, because of 
the small range of allowable impact parameters. The peak luminosity could be higher for 


(Frank fc Rees 1976; Lightman fc Shapiro 1977 
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tidal disruption events which occur in the orbital plane of a Kerr black hole. 

Following Rees (1988) and Evans and Kochanek (1989), the duration of a flare may be 
calculated based on the time that the infall rate is super-Eddington with the assumption 
that the circularization and accretion are rapid. The duration is 


f ^ f ■ 

^sup ~ ^mm 


/ JD \ —9/5 

18 r-®/- 1.5 


/ p \ 6/5 

1.9 r^/^m^'^^eg/^Mg’^^^years. 


(16) 

(17) 


If the viscosity is high enough that the matter has been accreted roughly as fast as it 
circularizes, the emission should fall off after tgup- 


5. Flare Energy Spectrum and Absolute Magnitudes 

The peak bolometric luminosity of the flares is determined from the scaling arguments 
given above. According to this scenario, the brightest flares are those in which accretion 
occurs in a thick disk (or in a slim disk), so a study of the spectrum for that case alone is 
the main focus. Cannizzo, Lee, and Goodman (1990) have investigated the spectrum a flare 
if the debris forms a thin disk. At the earliest and brightest times, the spectrum was found 
to peak at ~ 300A, with U and V absolute magnitudes of -16.7 and -15.2 for a IO^Mq black 
hole. 

Consider a thick disk spectra to begin with. Figure || shows the spectra calculated 
using a thick disk code (Ulmer, 1997). The spectra presented are typical thick disk spectra 
for a solar mass star disrupted at ~ 2Rt, the disruption radius (Eq. 0, by a IO^Mq black 
hole. Viewing angle will have a large effect on the spectra, because when seen edge-on, the 
bright, hot funnel is obscured by the disk. Viewed face-on, the disk is very bright, radiating 
slightly above the Eddington limit. When viewed edge on, the bolometric luminosity is less 
by a factor of ~ 50. For the thick disks considered here, the luminosity is 2 to 3 times the 
Eddington limit. The absolute magnitudes for a disk around a IO^Mq black hole radiating 
at 2.5 times the Eddington limit, are for face-on and edge-on disks: U:-21, V:-19.5 and 
U:-19, V:-17.5, respectively. For a similarly bright disk around a 10 ®Mq black hole, the 
absolute magnitudes are U:-17, V:-15.5 and U;-15, V:-13.5. The bolometric corrections are 
~ —7.5 mag in V, where reddening, which may well be important in the center of a galaxy, 
and dimming from cosmological effects (though in a flat cosmology, very high redshift flares 
would brighten as (1 + 2 :) due to K-corrections) are neglected. The disks may radiate closer 
to the Eddington limit than do the disk models above which radiate at 2 to 3 times the 
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Eddington limit. 


A thick disk may drive a strong wind. If a wind were dense enough, the effective 
photosphere would be moved outwards from the surface of the disk, and the spectra would 
appear redder. Some variable stars, namely luminous blue variables, show what are believed 
to be wind induced color changes, resulting in V-band brightening of 2 magnitudes (e.g. 
[Humphreys fc Davidson 1994|) . In super-Eddington thick disks, strong winds are to be 
expected on the purely phenomenological grounds that in stars, winds increase as a star’s 
luminosity approaches the Eddington limit. Additionally, instabilities in the funnel of thick 
disks may drive strong winds as well as jets (e.g. jJaroszyhsky, Abramowicz, fc Paczyhski 


198U|; [Mityananda fc Marayan 1982|). Alternatively, a strong wind could be driven before the 


debris completely settles into a disk. Rees (1988) observed that the total energy required 
to eject the majority of the bound material is small (see Eq. ^ ~ MqC^R^,Rs/R^. 


Limits on the strength of the wind are now examined. The mass loss rate is limited by 
the energy required to lift the material from the disk and carry it to inhnity. In the context 
neutron stars, Paczyhski & Proszynski (1986), have considered energy balance limits on 
winds. Following Paczyhski & Proszynski, the mass loss rate is 


m 




(Lb — -hph)2Rb 

TEdd2Rb 

c^s 


(18) 

(19) 


where Rh is the base radius from which the matter is lifted into the wind by the radiation, 
and Lb and Lph are the luminosity at the base radius and at the photosphere. If the base 
radius is taken to be ~ 10i?s, m ^ O.OSMQMgyears”^. The wind may, therefore, carry away 
a sizable fraction of the disk mass for a IO^Mq black hole. 


Given a maximal mass loss rate, the maximum photospheric radius may be estimated 
as follows. Near the photosphere, m = dyrnpRpb and pRph ~ rrip/aT', where rrip is the proton 
mass, and ctt is the Thompson cross-section. If n ~ Ugsc ~ ^JGMbh/R ph, then 

^ rh^cxT 

^ 327r2GMBHm2 

^ 120 Rs(—) (21) 

\VescJ 

The minimum effective temperature is expected to be ~ 1.1 x 10^ K, which is still very 
hot. The velocity of the outgoing wind could be lower than the escape velocity, and if the 
luminosity is very close to the Eddington limit at the photosphere, the escape velocity would 
be reduced by 1 — L/L-Edd from the expression given above, creating a larger photosphere. 
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Including a photosphere with T ~ 10^ K which radiates at the Eddington-limit in the thick 


disk spectra calculated above, damps the high energy spectrum but has little effect on the 
U and V band luminosities. The optical bands are weakly affected even though a larger 
fraction of the light goes into optical, primarily because the total luminosity of the source is 
diminished as energy is carried away by the wind and because a significant fraction of the 
optical light comes from part of the disk outside the photosphere at a few hundred R^- The 
optical luminosities of a thick disk and of a thick disk with a wind are nearly the same. 
The main effects of a strong wind are the reduction of the very high energy spectrum of the 
flare and a possible shortening of the flare by loss of the stellar debris to the wind. Such 
strong winds would greatly reduce the expected number of soft X-ray tidal disruption flares 
(Sembay & West 1993) accessible, for instance, to ROSAT. 

A third possibility, and the most optically luminous one, is that an extended, static 
envelope could be formed around the disk ( |Loeb &: Ulmer 19971) . Such an envelope might 
form either if the initial collisions of the debris streams or subsequent evolution of the disk 
stream expels much of the matter to large radii, and strong radiation pressure disperses 
this marginally bound gas into a quasi-spherical configuration. In this scenario the object 
is spherically symmetric with Eddington luminosity and effective temperature 



( 22 ) 


where Mgnv is the mass in the envelope. In this case, the bolometric corrections would be 
much less (~ 2.5), and the objects would be as bright as Seyferts. 

In addition to observing the very blue continuum of the thick disk, it may be possible 
to observe a flare in emission-lines. The disk may radiate in emission lines, but to what 
extent is uncertain because the emission would depend on the temperature and density 
structure both of the disk photosphere and of any wind. As discussed in Roos (1992) 
and Kochanek (1994), prior to circularization, the bound material generally has densities 
higher than lO^^atoms cm“^, so most of the energy the material radiates (or absorbs and 
re-radiates) will be in the continuum and not in emission lines. Moreover, such debris 
has a relatively small covering fraction. Recombination lines, e.g. Hq, might be radiated 
from the recombination of the debris or from the ionization and subsequent recombination 
of a neutral surrounding environment. The possibility that the disk may radiate in 
recombination lines was suggested by Eracleous et ah (1995) and Storchi-Bergmann et ah 
(1995), to explain observations of variable Ha lines in spectrum of the Seyfert NGC 1097. 
However, the disk is probably extremely optically thick even at thousands of Rg, with 
typical surface density O.5M0/(7r(lOOOi?s)^) of 3000Mg“^gcm“^, so any emission lines would 
have to be formed close to the surface of the disk and would be suppressed from optically 
thin estimates (Storchi-Bergmann et ah (1995), Eq. 2). 





6. Long-term Emission 


Because much of the flare’s energy is radiated between 10 and 100 eV, the possibility 
that the HII region created by the ionizing radiation of a flare could be observed is intriguing 
(e.g. Rees 1997). Assuming that ~ 5% of the rest mass of the star is radiated in the flare 
and using the thick disk spectra discussed above, the total number of photons between 
~ lOOA and lOOOA is ~ 10®^. The ISM surrounding the BH will be in a steady state, 
because the recombination time, tree ~ 10 ®n“^T 4 ^^years, is longer than the time between 
disruptions, (tdis ~ lO^years, where n is the number of atoms per cm^. The ionizing photon 
rate is therefore ~ 3 x 10^^s“^ or ~ 1% the total ionizing photon rate from the Galaxy 
(Mezger 1978; Bennett et ah 1994). Direct detection of from tidal disruption ionization 
is therefore unlikely. 


At higher photon energies, the tidal disruption events may (in a time-averaged sense) 
contribute the majority of photons. The production rate of He II ionizing photons by tidal 
disruption events is ~ 3 x 10®°s“^, which if a large fraction (e.g. 10%) of the photons ionize 
He H and ~ 20% of these ionization produce He H 4686, then the luminosity in He H 4686 
(Pacschen-line) would be ~ 2 x lO^^ergs s“^. The number of He H photons would be 
reduced by a factor of ~ 50 in the case of the strongest winds described above. For densities 
of one atom per cm^, the radius of the nebula would be 

1/3 


Rffli 


/ 3N^ 


1.5 X lO^^n ^/^T4'^^cm, 


(23) 


VdTrn^aB 

where cib is the case B recombination coefficient for hydrogen. A prediction of tidal 
disruptions is therefore a highly-ionized nebular emission in the inner kpc of galaxies. Such 
nebular emission may be produced by a variety of sources. He H 4686 lines for instance 
may be produced by Of stars ([Bergeron 1977|) or Wolf-Rayet stars (e.g. jSchmutz, Leither^ 
&; Gruenwald 199"^ , although such lines would generally be broad (~ 3 — lOA) whereas 
nebular lines would be narrow. Possible competing ionizing sources include shocks ( ^hull 
jfc McKeej) , X-ray binaries ([Garnett et ah 19'90D , and Wolf-Rayet stars. While all of these 
objects are concentrated towards galactic centers, it may be possible to distinguish between 
a tidal disruption nebula and most other ionizing sources. The fossil nebulae may be most 
readily observed in galaxies where the background ionizing level is small, corresponding to 
a small star formation rate, but since ample gas is required to absorb the radiation, the 
best candidates for detection of the nebulae may be Sa galaxies. Such nebulae may also be 
visible in radio lines. 


The idea that the unbound debris may produce a supernova remnant has been discussed 
by Khokhlov and Melia (1996). For typical disruptions, the unbound material is emitted 
in a fan of small opening angle, G < 0.1 rad^, with an energy of ~ 10^° ergs. Rare events 
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with pericenters close to the Schwarzschild radii, may provide 10^^ ergs. Such events may 
produce supernova remnants that could be observed in external galaxies. For remnants 
which sweep up 100 times their initial mass, the remnants would be displaced from the 
black hole by 

i?disp ~ (24) 

A supernova remnant lasts ~ lO'^ years, corresponding to one such remnant per galaxy. 


7. Conclusions 

Tidal disruption flares occur in two classes - those which are sub-Eddington and those 
which radiate near the Eddington limit. Flares from black holes above ~ 2 X IO^Mq will 
generally not radiate above the Eddington limit. 

The brightest flares will occur around ~ 5 x IO^Mq black holes, but they will be 
relatively rare as they require a hnely tuned impact parameter. The majority of the 
brightest events will occur around ~ 2 x IO^Mq black holes. The corresponding timescale 
for decay of such flares (Eq. is about one month. For all of the spectral scenarios 
considered, the optical luminosity is high with peak apparent magnitudes in U and V of 
-21 and -19.5 for a IO^Mq black hole and -19 and -17.5 for a 10 ®Mq black hole, neglecting 
reddening. Both high redshift supernovae searches and low redshift supernovae searches 
may observe the flares. As tidal disruption events occur at a rate of ~ 100 times less 
frequently than supernovae, many supernovae must be observed before a flare is expected 
to be seen. Furthermore, if the flares come from the centers of galaxies, then they may be 
harder to detect against the bright nuclei. The extreme blueness of the flares may help 
separate them from the background. In fact, because the spectra are so blue, K-corrections 
may actually brighten the flares in optical bands. 

The He II ionizing radiating produced in the flares may dominate that which is 
produced in the centers of quiescent galaxies, creating a steady state, highly ionized diffuse 
nebula with extent of ~ 1 kpc which may be observable in recombination lines in all galaxies 
with both tidal disruption events and gas near the nucleus. 
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Fig. 1.— Schematic of luminosity versus mass accretion rate (where MEdd = -hEdd/ec^), 
where e = 0.1 is the efficiency of the accretion process. The brightest accretion disks 
have high mass accretion rates and are thick. Thick disks are not much brighter than 
the Eddington luminosity, because at very high accretion rates pressure gradients push the 
inner edge of the disk closer to the marginally bound orbit, and energy is advected (because 
of geometry) (e.g. Paczynski and Wiita, 1980). Thin disks, (e.g. Shakura & Sunyaev 1973) 
and optically thin advection dominated models (e.g. Narayan & Yi 1995) are less luminous. 
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Fig. 2.— Timescales as a function of black hole mass a solar-type star with a pericenter 
equal to the tidal radius. For a IO^Mq, black hole, thick disks can form and bright flares can 
occur if face is less than ten which requires only that a > 1.5 x 10“®. Above ~ lO^M©, the 
Eddington time becomes shorter than the circularization time, so the disk cannot be thick, 
and the event will be less luminous. 






Fig. 3.— The fraction of tidal disruption events with peak accretion rates greater than 
the Eddington-limit are shown (upper curve) for the case of solar-type stars with a uniform 
distribution of pericenters and (lower curve) stars with a Salpeter mass function between 
O.SMq and IMq, also with a uniform distribution of pericenters. 
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Fig. 4.— Predicted spectra of a tidal disruption flare around a IO^Mq black hole. The flux 
is appropriate for a source at a distance of 100 Mpc. The two solid lines show calculated 
spectra from an optically and geometrically thick disk as viewed from face-on and edge-on. 
The former is much brighter and hotter because a large fraction of the energy is radiated in 
the funnel. The dashed line is a Planck curve with T = 10^ for comparison. 




